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In this paper, we investigate the single mode quantum properties of the codi- 
rectional Kerr nonlinear coupler when the frequency mismatch is involved and a 



condition for an exact solution of equations of motion is fulfilled. Particularly, we 
investigate quadrature and principal squeezing, Wigner function, quadrature dis- 
tribution, phase distribution and phase variance. We show that the quadrature 
squeezing and the phase variance can exhibit collapse-revival and collapse-revival- 
subrevival phenomena, respectively, based on the values of the detuning parameter. 
Furthermore, we analytically demonstrate that the system can generate cat states, 
in particular, Yurke-Stoler states. 
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I. INTRODUCTION 

Recently, there has been a great interest in the possibility of using optical devices for ultra- 
ligh-speed data processing. This is the main object in the quantum information theory, which 
aims at storing and transferring data Jj. One of the promising devices for data transmission 
is the nonlinear directional coupler that consists of two or more parallel optical waveguides 
fabricated from some nonlinear material. Both waveguides are placed close enough to permit 
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flux-dependent transfer of energy between them by means of evanescent waves. This flux 
transfer can be controlled by the device design and the intensity of the input flux as well. The 
outgoing fields from the coupler can be examined as single or compound modes by means 
of homodyne detection to observe squeezing of vacuum fluctuations, or by means of a set 
of photodetectors to measure photon correlations, photon antibunching and sub-Poissonian 
photon statistics in the standard ways. The investigation of the quantum properties of 
light propagating in the directional couplers have attracted much attention for generating 
nonclassical light (see the review paper |2J and the references therein). Moreover, directional 
coupler is experimentally imp 
and certain organic polymers 

Among the different types of directional couplers the directional Kerr nonlinear coupler 
(DKNC) has taken much attention as a result of its application in optics as an intensity- 
dependent routing switch 0, Q]. DKNC is useful for low intensity fields and periodical 
exchange of energy between the guides; but for high intensity fields, energy is trapped by 
nonlinearity in the guide into which it was initially launched (self-trapping effect). T he q uan 



emented, e.g. in planar structures |3(, dual optical fibres 
3- 



tuin properties of the DKNC have boon studied by several authors '{. 3, 1Q, H, 12. llJ. fl^}]. 
For instance, in [lOj quadrature squeezing, principal squeezing and integrated intensity vari- 
ances have been calculated in an analytical way and investigated based on the transmission 
of light between waveguides. Moreover, quantum statistics for contradirectional KNC have 
been investigated [yj] showing that the switching between waveguides is accompanied by 
nonclassical effects, e.g. two-mode squeezing, the generation of pure state and single mode 
photon antibunching. Influence of the geometry of the coupler (when the linear coupling 



coefficient is variable) on the statistics of the DKNC ll|, |l5j (also contradirectional KNC 



|15l |) was investigated showing that there is a possibility to control energy switching between 
waveguides by adjusting the form of a coupling function. The numerical technique based 
on the diagonalization of the Hamiltonian is used for studying the mean-photon number 
[i| and the phase distribution in the framework of a quasiprobability distribution func- 
tion Jl^^Furtherniore, in the most of papers dealing with DKNC in the literature, e.g. 

[13I . [14J], authors have shown that the evolution of the mean-photon numbers 
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exhibit collapse-revival phenomenon arising from the nonlinear exchange of energy between 
waveguides. In this regard there is a similarity between DKNC and the behaviour of the 
atomic inversion in the Jaynes-Cummings model 16j. In the present paper we provide-under 
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a certain condition-the exact solutions for the equations of motion of DKNC taking into 
account the frequency mismatch. This would lead to many important results additionally 
to what are already known for DKNC. For instance, we show that the quadrature squeezing 
can exhibit collapse-revival phenomenon rather than the mean-photon numbers, as shown 

nnnnn 

earlier [7], [KJ, [Hi, ll3|, [14j . Also the phase variances exhibit collapse-revival-subrevival phe- 
nomenon similar to that of the two-mode Jaynes-Cummings model [17]. Moreover, for the 
first time, we analytically and numerically prove that the system can generate Yurke-Stoler 
states at certain times despite the coupling of the two modes. 

Now the Hamiltonian controlling the codirectional Kerr nonlinear coupler is given as 

H 2 

— = ^2[ujjOjaj + xafaj] + X^i^l^ + K(a\a 2 + a\ai), (1) 
i=i 

where uj\ and LO2 are the frequencies of the first and the second modes with the annihilation 
operators a± and a 2 , respectively, % an d X are the coupling constants proportional to the 
third-order susceptibility x an d responsible for the self-action and cross-action processes, 
respectively, k is the linear coupling constant between the waveguides. Kerr coupler can be 
implemented from certain organic polymers with high third-order nonlinearities Also 
cento-symmetric optical fibres can be adopted. We proceed that in above mentioned papers 



related to DKNC, i.e. 



10, 



11 



12 



13 



14j , the authors have neglected the non- relevant terms, 



i.e. the nonlinear rotational terms, for obtaining a closed form solution of the equations of 
motion related to (JT]). In this paper we restrict ourselves to the case in which the equations 
of motion can be solved exactly. In this case the system is linear in the sense that the mean- 
photon numbers exhibit oscillatory behaviour indicating periodic energy exchange between 
the waveguides even though there is a nonlinear medium between the waveguides. Actually, 
this does not mean that the coupler cannot generate nonclassical effects. More illustratively, 
the mean-photon numbers are given by the diagonal elements of the density matrix only. 
Nevertheless, the quantities, which depend on the off-diagonal elements, such as quadratures 
squeezing, Wigner function, phase distribution, etc., may generate nonclassical effects. Thus 
we are going to investigate the behaviour of these quantities in the present paper, considering 
also the influence of frequency mismatch A = uj\ — 002 ■ It is worth remembering that the 
phase mismatches were considered to the other types of coupler, e.g. linear coupler |l8| . 
nonlinear coupler [3] and Raman- Brillouin couplers [20I . 

We conclude this section by writing the equations of motion for ([1]) under the condition 
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Oi(*) = exp(-«At/2)|ai(0) [cos(At) -z^sin(At)] - if a 2 (0) sin(At) 
a 2 (£) = exp(-iAt/2)|a 2 (0) [cos(At) -z^sin(Ai)] - if oi(0) sin(At)} 



X = 2x, 

— — = —ioj\a\ — £%x\a\a\ + a^aij&i — IKCL2, 

Lib 

^-^ = — itu 2 a 2 — 2ix(a[a 1 + a 2 a 2 )a.2 — i/cai- (2) 

One can easily check that ajai + = C* is a constant of motion and hence the general 
solution for (ED is 



(3) 



where A = y k 2 + ^A 2 , A = u\ + 102 + 4x(7 and A is the frequency mismatch. One can 
easily check that d\(t) <-> 02(t) when ai(0) <-> 0.2(0). Thus we restrict the discussion to 
the behaviour of the first mode only. It is worth mentioning that for solving the problem 
in the space domain we have to use the substitution z = t$, where is the velocity of 
light in the waveguide and z is the travelled distance. Moreover, we do not consider the 
dissipation, which generally leads to decrease of the total number of photons and to the 
tendency to reduce the exchange of photons between the waveguides, i.e. to the reduction 
of the nonclassical effects @, 9]. 

Finally, using ([3]) we investigate only the single mode quantum properties for the DKNC 
We perform such investigation in the following order: In section 2 we examine the quadra- 
tures and principal squeezing. In section 3 we calculate and discuss the quasiprobability 
distribution and quadrature distribution. In section 4 we study the phase distribution and 
its variance. The results are summarized in section 5. 



II. QUADRATURES AND PRINCIPAL SQUEEZING 

The photons produced in a nonlinear optical device such as DKNC are known to have 
unusual correlation properties, which results in many nonclassical aspects of the radiation 
field. Thus in this section we demonstrate the single mode quadratures and principal squeez- 
ing. As is well known squeezed light has less noise than coherent light in one of the field 
quadratures. This light can be measured by homodyne detection where the signal is super- 
imposed on a strong coherent beam of the local oscillator. Moreover, squeezing is one of the 



5 



most important phenomenon in quantum optics because of its applications in various areas, 



2l| 



e.g., in optics communication, quantum information theory, etc. 

To investigate the single mode squeezing we define two quadratures X and Y, which 
denote the real (electric) and imaginary (magnetic) parts, respectively, of the radiation 
field, as 

X A{{t)], Y = -^=[ix(t) - A\(t)}, (4) 

where Ax(t) = &i(f) exp[^(u;i +0^)]. As we mentioned in section 1 we restrict the discussion 
to the first mode. The quadrature operators (jl]) satisfy the following commutation relation 



X,Y 



Therefore, the uncertainty relation related to (jSJ is 



((AX) 2 )((AF) 2 )>-, 



(5) 



(6) 



where the variance ((AX) 2 ) = (X 2 ) — (X) 2 and similarly form can be given for ((Ay) 2 ). 
We say that the system exhibits X-quadrature squeezing when 

S = 2((AX(f)) 2 ) - 1 < 0. (7) 

The equality sign holds, i.e. S = 0, for minimum-uncertainty states. Similar definition can 
be given for the F-quadrature (defining a Q-factor). 

On the other hand, the quadratures (jl]) can be represented by the Hermitian operator 



V2 



[A x (t) exp(-i0) + A\(t) 



(8) 



where is a phase, which can be controlled by the homodyne detector. The operator (JHJ) 
reduces to the X-quadrature (the in-phase component of the field) when = 0, while when 
(j) — 7r/2, it gives the F-quadrature (out-of-phase component of the field). According to (IE|) 
the variance is dependent on t he p hase 0, which can be controlled to give the minimum 



of a 



22 



1 quadrature variances 



22 



23| . Thus we obtain the notion of the principal squeezing 



231 ]. which can be expressed for the single-mode as 

i#) = 2[(ii(f)i 1 (f)) - (AimAM) - \(Ai( t )} - (A&n 



(9) 



Squeezing occurs when 77(f) < 0. 
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FIG. 1: Evolution of the squeezing factors and principal squeezing of the first mode when k = 
l,X = 0.5, («i,a 2 ) = (0.2,0.2) and for (a) A = 0s" 1 (solid curve for S(t), short-dashed curve for 
Q(t) and long-dashed curve for rj(t)), (b) evolution of the factor Q(t) for A = 50s . The straight 
line in (a) is given to show the squeezing bound. 



The different moments of the operator A\(t) can be evaluated as 



{A\ m (t)A^t)) = a?(t)ar(*)* [ * (n-1) ~^ (m_1)] 



exp[£[z 1 



1)], 



where 



and 



|«i| 2 + |a 2 | 2 , z = exp(—2ixt), ai(t) = a x (t) + ia y (t) 

r \ t \ \ i \ r ^ ,sin(At) 

a x {t) = oti cos(At), a y (t) = —[ai— + a 2 K,\ 



(10) 



(11) 



(12) 



2 ' X ' 

where a\ and a 2 are the initial field amplitudes. Throughout the paper we consider that 
«i and «2 are real. Expression (ITU]) reflects two facts: (i) When \ — the system reduces 
to the up conversion process, which switches the energy between the waveguides without 
generating nonclassical effects, (ii) When n = m all the moments are independent of x and 
the mean-photon numbers exhibit oscillatory behaviour only rather than collapse-revival 
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11 



12, QQ. Moreover, this fact 



pattern, which is representative for the DKNC 
shows that the system exhibits always Poissonian statistics. 

Now we start the discussion by investigating the behaviour of the quadratures squeezing. 
We show that the system can provide squeezing and the quadrature squeezing can exhibit the 
collapse-revival phenomenon. Furthermore, squeezing cannot be simultaneously generated 
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in the two quadratures. In doing so we use (flUI) to express the squeezing factors S(t) and 
Q{t) as 

S(t) = 2\a 1 (t)\ 2 + G 1 (t)-G 2 (t), 

(13) 

Q(t) = 2\a 1 (t)\ 2 -G 1 (t)-G 3 (t), 

where 

G x {t) = 2[(«2( t ) _ cos(e(2 X i)) + 2« x (t)^(t) sin(e(2 X t))]/(2 X t), 

G 2 (t) = 4[a x {t) cos(e(xt)) + a„(t) sin(e( X i))] 2 / 2 (xt), 

G 3 (t) = 4[^(t) sin(e( X t)) - a„(t) cos(e( X t))] 2 / 2 (xt), ( 14 ) 
e(nxt)) = n(n — l)%t + e sin(2nx^), 

/(nx*) = exp[-2esin 2 (nxt)]- 
It is obvious that the behaviours of S(t) and Q(t) are periodic as a result of the nature of 
the coupler, which basically depends on switching energy between waveguides by means of 
the evanescent waves. Moreover, the origin of the occurrence of the nonclassical effects is 
in the third order nonlinearity, which is related to the envelope function f(nxt), i.e. to the 
nonlinear phase modulation term, that causes such effects. Furthermore, involving squeezing 
factors, the function f(nxt) indicates that they can exhibit collapse-revival phenomenon, 
as we shall see. Furthermore, the generation of squeezed light in DKNC is quite obvious 
from (CQ), which-apart from the free part-includes two main parts: linear- interaction part 
and self-cross nonlinear interaction part. The latter is well known in the literature, see e.g. 



24j . as being able to generate nonclassical light such as squeezed light, whereas the former 
switches energy only. Now we discuss some analytical results for the system by focusing the 
attention on a simple case, when A = 0s _1 , a,\ = a 2 — a. In this case expressions ( |T3l) can 
be written in the forms: 

S{t) = 2a 2 + 2a 2 cos[2tA + e(2 X t)]/(2 X t) - 4a 2 cos 2 [iA + e( X t)]/ 2 ( X t)], 

(15) 

Qit) = 2a 2 - 2a 2 cos[2tA + e(2 X t)}f(2 X t) - 4« 2 sin 2 [tA + e( X t)}f 2 (xt)]- 
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FIG. 2: Evolution of the squeezing factors and principal squeezing of the first mode when k = 
l,X = 0.5, (ax,a 2 ) = (2,0) and for (a) A = Os" 1 (solid curve for S(t), short-dashed curve for Q(t) 
and long-dashed curve for r](t)), and A = 50s" 1 for Q(t) (b) and r/(t) (c). The inset in (a) shows 
that Q(t) cannot exhibit squeezing, while that in (b) shows the amount of squeezing, which can be 
obtained. 

It is easy to prove that 

S(t) + Q(t)=Aa 2 [l-f 2 ( X t)}. (16) 

It is evident that expressions (jT5l) provide extreme values when the envelope functions tend 
to unity. This means that when the system generates squeezed light S(t) + Q(t) ~ and 
then squeezing cannot be simultaneously generated in the both the quadratures. 
The envelope function in the second term of ( |T5l) is maximum when 

x t = m'ir/2, m' = 1,2, ... (17) 

In this case expressions (fT5l) reduce to 

S(t) = 2a 2 {[l - (-1)'"'] - 2[/ 2 (m'f ) - (-1)™'] cos 2 (tA)}, 

(18) 

Q(t) = 2a 2 {[l - (-1)'"'] - 2[/ 2 (m'f ) - (-1)™'] sin 2 (tA)}. 

From ([TBI , for these particular values of the interaction time, it is easy to prove that 
squeezing occurs for weak-intensity regime, i.e. < a < 1, where /(m'|) ~ 1 for m! = 
1, 3, 5, • • • . Nevertheless, squeezing cannot occur for strong-intensity regime, i.e. when a > 1 
where f(xt) — 0. On the other hand, the third term in (fT5l) is maximum when 



Xt = m'ff, m' = 1, 2, 



(19) 
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FIG. 3: The single mode W function for k = l,x = 0.5, A = 0s 1 ,t = tt and (a) (a>i,a2,D) = 
(2,0,0), (b) (2,2,4). 



(b) 




FIG. 4: The single mode W function for k = l,x = 0.5 and (a) (A, ai, at2, t) = 
(0s _1 , 0.2, 0.2, 3.139997), (b) (0s" 1 , 2, 0, 6.36005), (c) (50s" 1 , 2, 0, 6.36005). 

In this case ([TBI) reduces to S(t) = Q(t) = 0, i.e. the system generate minimum-uncertainty 
states. In spite of this fact we have numerically found that the system can generate squeezing 
close to xt = m 'K regardless of the values of atj. This is related to the factor e(nxt) involving 
in the trigonometric functions. 

All these analytical facts and the influence of A on the evolution of quadrature squeezing 
are presented in figures 1 and 2 for given values of the parameters. From Fig. 1(a), where 
the intensities are weak, one can see that squeezing periodically (but not simultaneously) 
occurs in both the quadratures. The influence of the detuning parameter is given in Fig. 
1(b), which displays an exact periodic collapse-revival phenomenon. This is related to the 
values of the scaled interaction time tX, which is rapidly changed when A is large. More 
illustratively, the squeezing factors include two forms of periodic functions: one is coming 
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from the self-cross- nonlinear interaction part, in particular, the envelope function whose 
period is n/x, and the other is arising from the linear- interaction part whose period is n/X. 
Thus when the values of A increase, the period of the energy exchange between waveguides 
decreases, i.e. many oscillations occur, till the interaction time becomes t = tt/x at this 
moment the field is trapped instantaneously by nonlinearity in the waveguides and the 
squeezing factors show collapse. As the interaction proceeds the phenomenon is periodically 
repeated. Now we draw the attention to Figs. 2. From Fig. 2(a) it is obvious that the system 
produces vacuum and coherent light periodically. This can be realized analytically since for 
A = 0s _1 ,t = (m' + l/2)7r, A = 1 the amplitude |cm(i)| = (regardless of the values of the x) 
and the system produces the vacuum state. Similar arguments can be given for the coherent 
light. Also from Fig. 2(a) squeezing occurs in the S(t) only (see the solid curve in the inset) 
at particular values of the interaction time. In fact, at these values the system generates 
Yurke-Stoler coherent states (YSCS) 25], as we shall show in sections 3 and 4. Moreover, 
squeezing can be established in Q(t) when the frequency mismatch is included (see Fig. 
2(b)). Also from Fig. 2(b) Q(t) exhibits particular type of collapse-revival phenomenon, 
which can be explained in the following sense. According to the values of the interaction 
parameters considered in Fig. 2(b) (i.e. a% = 0, A ^> Is -1 and A ~ A/2) expressions 
( TT5|) — (TT9T) and the discussion around can be used. In this regard Gi(t) exhibits two times 
revival patterns compared to those occurring in G^(t). When revivals occur in G\(i) and 
Gz(t) simultaneously, they destructively interfere and cancel out each others showing such 
shape, however, at this stage the system generates squeezing, i.e. squeezing occurs close to 
X = mV. The comparison between figures 1(b) and 2(b) shows that as the values of otj 
increase the widths (i.e. the envelopes of the revival patterns) decrease, but the collapse 
period is enlarged. Actually, we found that S(t) provides quite similar behaviour as Q(t). 

Now we turn the attention to principal squeezing i](t), which for the mode under consid- 
eration can be expressed as 



T)(t) = 2|a 1 (t)| 2 |l - exp[-4esin 2 (xt)] - exp[-4£ sin 2 (xt)] 

x [l + exp[-8e sin 2 (xt) cos(2xt)] - 2 exp[-4e sin 2 (xt) cos(2xt)] cos[2xt - 4e sin 2 (xt) sm(2x<)]] 
As we did for the quadratures squeezing one can prove that for jt = mir the principal 
squeezing vanish, however, for yt = m'ir/2, m' is odd integer, it reduces to 
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r } (t) = -4\a 1 (t)\ 2 exp(-4e). (21) 

This means that squeezing always occurs provided that atj are finite. This is related to 
the fact that the phase of the homodyne detector is adjusted for obtaining the minimum 
variance. For the resonance case principal squeezing factor r)(t) is presented by the long- 
dashed curves in Fig. 1(a) and Fig. 2(a) for weak and strong intensities, respectively. 
From Fig. 1(a) it is obvious that rj(t) provides the envelope for the periodic nonclassical 
effects in the Q(t) and S(t) factors. More precisely, rj(t) starts from zero before switching 
on the interaction, monotonically decreases (increasing squeezing) as the interaction time 
increases till yt — 7r /2 showing its minimum (maximum squeezing) and hence increases 
monotonically till providing its initial value at ■yt = 71 ■ This behaviour is periodi cally 
repeated. Similar behaviour has been seen for the anharmonic oscillator model 22j, |23j. 
From Fig. 2(a) for strong intensity regime rj(t) exhibits behaviour, which is similar to that 
of S(t). Nevertheless, rj(t) provides nonclassical squeezing over intervals of the interaction 
time larger than those for S(t) (compare the long-dashed and solid curves in the inset). 
Now we draw the attention to the non- resonance case, which is presented by Fig. 2(c) for 
given values of the interaction parameters. We have noted that r](t) cannot provide collapse- 
revival phenomenon. Moreover, for weak intensity regime rj(t) is insensitive to the value of 
A, which almost leads to the long-dashed curve in Fig. 1(a). From Fig. 2(c), which is given 
for strong-intensity regime, we can see that the behaviour is completely different from that 
of the Q(t), i.e. it is periodic involving nonclassical effects (squeezing) rather than collapse- 
revival phenomenon. Surprisingly, Fig. 2(c) is similar to Fig. 2(e) 10( for the resonance 
case. 

Generally, we can conclude that the origin of the occurrence squeezing and collapse-revival 
phenomenon in the single mode quadratures squeezing lies in the competition between the 
nonlinearity and the frequency mismatch. The locations of the revival patterns in the time 
domain depend on the values of t\ and their shapes depend on the intensities of the field 
launched in the waveguides initially. Finally, the behaviours of the quadrature squeezing 
presented here are completely different from those given in [h]] as a result of neglecting 
the nonrelevant terms there. Also squeezing in the framework of the principal squeezing 
has been remarked. Nevertheless, the principal squeezing cannot exhibit collapse-revival 
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3.50 



FIG. 5: The quadrature distribution P(x) of the first mode for k = 1, x = 0.5 and for (a) t = ir,cti = 
2 and (a 2 ,A 2 ,D) = (0,0s~ 2 ,0) (solid curve), (2,0s -2 , 4) (short-dashed curve) and (0,5s~ 2 ,0) 
(long-dashed curve); (b) (A,ax,a 2 ,t) = (0s~\ 0.2, 0.2, 3.139997) (solid curve), (0s~\ 2, 0, 6.36005) 
(short-dashed curve) and (50s -1 , 2, 0, 6.36005) (long-dashed curve). 

phenomenon. 



III. QUASIPROB ABILITY DISTRIBUTION FUNCTION 

Quasiprobability distribution functions (PF-Wigner, Q-Husimi and P-Glauber functions 
261 ]) are important tools to give insight in the statistical description of a quantum mechan- 



ical system. These functions can be measured via homodyne tomography [27j. Here we 
investigate the single mode quasiprobability distribution functions and quadrature distri- 
bution for the system under consideration using the technique given in [28[]. Furthermore, 
in the following section we use these functions in investigating the phase distribution and 
phase variance. In order to obtain these functions we have to calculate the single mode 
s-parameterized characteristic function having the form 

C((,t,s) = Tr {p(0) expfcito - CMt) + ^IC| 2 ]} , (22) 

where p(0) is the density matrix of the system and s is a parameter that takes on the values 
1, and —1 corresponding to normally, symmetrically and antinormally ordered characteris- 
tic functions, respectively. From (fTUj) and (|22j) the characteristic function can be expressed 
as 
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rti,n 2 =0 



x ^[^(„ 2 -l)-^( ni -l)] exp ^^„ 2 - m _ ^ 



(23) 



where we have considered that the two modes are initially prepared in coherent states with 
amplitudes a%, a%. 

The s-parameterized quasiprobability functions can be obtained through the relation 



W((3, M) = ^ / C(C> *, a) exp(/3C - P*0d 2 (. 
On substituting (l23j) into (124")) and using the technique given in 



(24) 



we arrive at 



W{/3,t,s) 



^- T exp( 2 1 "^ 



ni,n2=0 



(_i)"i a^ffla*" 1 (t) ( 2 \"2 
n 2 ! 



(25) 



X2 .[^(n 2 -l)-^(m-l)] e 3 t p[ e ( 2! na-n 1 _ l)]L^-»i 



where L™(.) is the associated Laguerre polynomial of order n and (3 = x + iy = \/3\ exp(iO), 
being the phase of (3. One can easily verify that when yt — and m' is integer, the 
system produces coherent light with amplitude cti(t), as we mentioned in section 2. This 
can be verified from (|25|) by means of the generating function for the Leguerre polynomials 
(see ( 147!) in the Appendix). In this case ( |25l) reduces to the W function for the coherent 
state 



W(j3,t,s) 



exp 



1-s 



|/5-«i(t)p 



(26) 



7T(1 - S) 

Furthermore, ( |26l) indicates that the distribution cannot provide negative values at the phase 
space origin, i.e. (3 — 0. Generally, (125]) includes complicated quadratic phase factor, which 
plays an essential role in generating cat states 29] and it reflects the strong entanglement 
between the two modes. In fact, for particular values of the parameter t\ expression (T25]) 
can reduce to that for the cat states. For instance, when t\ = (m' + 1/2)tt and ml is integer, 
( I25p can be modified to the following form: 



(27) 



W((3, t ,8) = ^\ exp - ^(t)! 2 ) + exp (-^-J/3 + ^(t)| 2 ) 



+ 



2 exp [-^(|/3| 2 + D)] sin (^(/Jaftt) + /3*a a (*)))}, 
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FIG. 6: The single mode phase distribution P(Q,t,s = —1) for k = l,x = 0.5, (A, ai, 012) = 
(0s -1 ,2, 2) and for t = (star-centered curve), 7r/4 (short-dashed curve), 2.94 (long-dashed curve) 
and 7r (solid curve). 

where 

D = e- |«i(t)| 2 . (28) 

The derivation of ff27j) is given in the Appendix. Moreover, one can easily verify that (|27|) 
is normalized, 

r W(P,t,s)d 2 p = I. (29) 



The form (127)) reduces to that of YSCS 25| when D ~ 0. In other words, when D ~ 
the system generates the states 



M*)> = +exp(z|)| -ia x (t))}, 



(30) 



which is YSCS. As is well known YSCS exhibit Poissonian statistics and provide squeezing 



30] . These results are obtained in section 2. This indicates that the cat states generated in 



DKNC are YSCS. Moreover, when A = 0s 1 and from (fT21) . expression (1281) takes the form 

D = «iSin 2 (At) + «2Cos 2 (At). (31) 

Expression (13T1) provides several consequences: (i) DKNC generates YSCS when one of 
the modes is initially prepared in the vacuum state and the other is in the coherent state 
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provided that the length of the waveguides (or the interaction time) is controlled so that 
D — 0. (ii) when A = 0, i.e. the linear interaction between the waveguides is neglected, the 
system can generate YSCS by choosing cti ^ and a<i = 0. (iii) the system cannot generate 
YSCS when «i = a 2 = even though D = 0, because the expression ( 1271) reduces to that 
for the vacuum. 

On the other hand, when D > 1 the interference part in (l2Tj) is destroyed and hence the 
DKNC generates the statistical-mixture coherent states having the form 

P(t) = 2 [l*ai(*))<i«i(*)l + I " iai(*))(-iai(*)|] • ( 32 ) 

An illustration of (|2"T|) (or fl25|) ) is given in Figs. 3(a) and (b) for D = and 4, respectively. 
In Fig. 3(a) the W function consists of two Gaussian bells, corresponding to the statistical 
mixture of individual composite states (cf. ( 130]) ) and interference fringes inbetween arising 
as a result of the contribution originating from the quantum superposition. In Fig. 3(b) 
the interference fringes are completely suppressed and the form of the statistical-mixture 
coherent states is well pronounced. Actually, by controlling the values of the A, the form 
of the cat states can be less or more pronounced. For instance, Fig. 3(a) can be reduced to 
Fig. 3(b) by taking A = ^/Es^ 1 instead of zero. 

Now we shed the light on the general case, e.g. apart from the case t\ = (m' + l/2)7r. 
This will be numerically done by figures 4 for the W function obtained from ( 1251) . The values 
of the interaction parameters in Figs. 4 have been selected such that the system provides 
quadratures squeezing (see Figs. 1 and 2). In all these figures we can observe the nonclassical 
effects, such as negative values, stretching, multipeak structure and deformation. In Fig. 
4(a), where the initial intensities of the modes are weak, the W function exhibits contour- 
stretching as well as deformation around the phase space origin. The origin of such behaviour 
is in the generation of cat states in the microscopic regime | [28l |. where the contribution of the 
different components of the cat are located close to the phase space origin competing each 
others. When the values of the intensities and the value of the interaction time are increased, 
the shape of the cat states becomes more pronounced by involving multipeak structure 
(asymmetric peaks) (see Fig. 4(b)). This is in a good agreement with the information 
given above even though t\ ^ (m' + l/2)7r. We proceed by drawing the attention to the 
influence of the detuning parameter on the behaviour of the W function, which is given in 
Fig. 4(c). In this figure the W function exhibits almost stretched single peak structure 
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as well as negative values indicating that the nonclassical effects are still occurring. The 
transition from Fig. 4(b) to Fig. 4(c) through changing the values of the detuning is related 
to the fact that in the former the amplitude a\(t) is real, while in the latter it is complex 
(one can check this for the chosen values of the interaction parameters). This leads to that 
the oscillations arising from the nonlinear phase modulation are superimposed with those 
of the field amplitude ai(t) causing such behaviour in Fig. 4(c). We conclude this part by 
mentioning that a connection between the Fourier coefficients and the Q function for the 



case 2x = X using the diagonalization approach is given in 12] . 



We close this section by investigating the quadrature (or position) distribution P(x,t) 



which can be measured in the homodyne detector 
evaluated via the W function through the relation 



251 ] . The distribution P(x,t) can be 



oo 

P{x,t) = J W(x + iy,t,0)dy. (33) 



Substituting fl25l) into fl33l and after lengthy calculation, we arrive at 

y— oo min(ni,n 2 ) i (ni+ „„) _„ 2 M -* ni M 

P(x,t) = ^eM-^) E E { - 2)2 t:J^r ® 

" ni,n2=0 r=0 



xz [^(n2-i)-^(m-i)] exp [ £ (^ 2 -n 1 _ i)]H ni _ r (V2x)H n2 _ r (V2a; 



(34) 



where H m (.) is the Hermite polynomial of order m. It is more convenient to give the explicit 
form for P(x, t) for the case t\ = (m! + l/2)7r 

P ( x ' *) = yfc {exp[-2(x + a y {t)f] + exp[-2(x - a y (t)) 2 ] 

(35) 

+2exp[-2(x 2 + a 2 {t) + D)]sin(Axa x (t)f 



In Figs. 5(a) and (b) we have plotted P(x,t) corresponding to the W functions shown in 
Figs. 3 and 4, respectively. In general one can observe that P(x, t) provides nonclassical 
effects by including oscillatory behaviour. In Fig. 5(a) the solid curve shows the oscillatory 
behaviour related to YSCS, which is a direct consequence of the interference in phase space. 
Nevertheless, the short-dashed curve is a Gaussian bell having its maximum value at x = 0. 
This indicates that the system generates vacuum light. Actually, for this curve ct y (t) ~ 
and the interference part is completely suppressed since D = 4. This may contradict to the 
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result illustrated in Fig. 3(b) for the corresponding W function, in which the generation 
of the statistical-mixture coherent state is obvious. This confusion can be removed if the 
attention is drawn to the momentum distribution P(y,t), which exhibits, for this case, a 
two-peak structure (we have checked this fact). The influence of the A is presented by the 
long-dashed curve, which provides the two-peak structure indicating the generation of the 
statistical-mixture coherent state. On the other hand, we can note that there is almost 
consistence between the behaviours in Figs. 4 and Fig. 5b, where P(x,t) provides two 
asymmetric peaks for strong intensities regardless of the values of A and a single peak for 
weak intensities. We conclude this section by mentioning in {31] stated that "the concept 
of interference in phase space readily explains the similarity between the photon number 
distribution and the position distribution". This is not always correct. For instance, P(x,t) 
of the YSCS provides oscillatory behaviour (see the solid curve in Fig. 5(a)) even though 
the photon number distribution of YSCS is Poissonian, as we mentioned in section 2. 

IV. PHASE DISTRIBUTION AND ITS VARIANCE 

In this section we study the evolution of the phase distribution and its variance. The 
phase distribution associated with the first mode can be obtained from the quasidistribution 
( 125]) by integrating W((3, t, s) over the radial variable as 

POO 

P(G,t,s)= / \p\W(P,t,s)d\P\. (36) 
Jo 

On substituting (125]) into fl36|) and carrying out the integration, we arrive at 

00 "1 (-i)"i+™r(m+ " 2 ~" 1 +i)|ai(t)|"2+"i 



P(e,f, a ) = f 1 + 2 £ E M^p^i«u ;i <_ t ,_2_y 

\ ' ' / 2n \ ^— ' ^— ' (ni—m)\(n2—ni+my.mi VI— s/ 

k ri2>n\ m=0 



x cos[y? + {n\ — n 2 )Q] exp[— 2e sin 2 ((n 2 — ni)xt)} 
where 

det(t) = \ai(t)\exp[i(f)(t)], 

(p — (n 2 — ni)4>(t) + [n 2 (n 2 - 1) - n\{n\ - !)]%£ + e sin(2xt(^2 - n\)) 



(37) 



(38) 



and T(.) is the Gamma func 



tion. The phase distribution given by (1371) is just a special case 



of the formula (2.59) in 32j. As is well known the phase distribution obtained from the 
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FIG. 7: The single mode phase variance ((A<3?i(i)) 2 ) against the interaction time t for k = l,x = 
0.5 and (a) (A,ai,a 2 ) = (0s-\2,2) and (b) (50s- 1 , 2, 2) . 

quasiprobability distribution functions includes some difficulties. For instance, the singular- 
ity of the P function and the negative values involved in the W function for some quantum 
mechanical systems may be reflected in the corresponding phase distributions (see, e.g., 



331]). These difficulties can be avoided by using Q function, which is well defined and is 
always positive, in evaluating the phase distribution. Illustration of P(Q, t,s = —1) is given 
in Fig. 6 for given values of the parameters. From this figure one can see how the cat states 
can be generated in the system depending on the values of the interaction time (or the 
length of the waveguides). For instance, at t = 0, i.e. before switching on the interaction, 
P(0) exhibits single-peak structure around G = 0, which is representative for the coher- 
ent light (see the star-centered curve). As the interaction proceeds, i.e. energy exchange 
between the waveguides starts to play a role, the initial peak reduces to two wings around 
O = ±7r, which becomes straight line at t\ = tt/2, i.e. the phase distribution exhibits 
normal distribution (this case has not been included in the figure). It is worth remembering 
that the system generates vacuum state in this case. When the interaction time becomes 
close to 7T the form of the cat states starts to appear as a broader two-peak structure (see 
the long-dashed curve in Fig. 6). Eventually, at t = 7r the distribution provides two Gaus- 
sian peaks around O = ±7r/2 indicating the generation of the cat state, in particular, the 
statistical-mixture coherent state (cf. (|32|) ). Actually, the phase distribution is insensitive 
to the interference in phase space, i.e. the phase distribution of the even, odd, Yurke-Stoler 
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and the statistical- mixture coherent states are almost similar [34|. This leads to the fact 
that the phase distribution is insensitive to the values of A, which can change cat states 
to the statistical-mixture coherent states (see the discussion given in section 3 for the W 
function). 

The single mode phase variance is defined as 

((A* 1 (t)) 2 ) = (%(t))-(* 1 (t)) 2 . (39) 
The Zth moment of the phase distribution can be evaluated from (137j) through the relation 

($[(t)) = J Q l P(Q,t,s)dQ. (40) 

— 7T 

The evolution of the ((Al>i(t)) 2 ) is given in Figs. 7(a) and (b) for A = 0s" 1 and 50s _1 , 
respectively. In Fig. 7(a) the periodic behaviour is dominant indicating switching of energy 
between the waveguides. Also one can see that the system provides its initial stage periodi- 
cally. The detuning parameter A reflects itself in ((A&i(t)) 2 ) as collapse-revival-subrevival 
phenomenon (see Fig. 7(b)). As we mentioned in section 2, when the values of the A in- 
crease, the period of switching of energy between the waveguides decreases, causing such 
phenomenon. Furthermore, the main revivals occur around the values of the interaction 
times at which the system reduces to its initial form, however, the subrevivals occur when 
the system generate cat states. Actually, the density matrix of the superposition states 
has different components, each of them has its own collapse-revival pattern (when A is 
large), which interfere with each others producing subrevivals. Thus we can conclude that 
the occurrence of the collapse-revival-subrevival phenomenon is a direct consequence of the 
generation of different types of states in the system, i.e. vacuum, coherent, and cat states. 
It is worth mentioning that such behaviour has been observed for the single mode phase 



variance of the two-mode Jaynes-Cummings model 
here is more systemic. 



17] . however, the behaviour presented 



V. CONCLUSION 



In this paper we have discussed the single mode quantum properties for the codirec- 
tional nonlinear Kerr coupler, when the frequency mismatch is involved. The attention is 
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focused on the case 2x = X f° r which the solutions of the equations of motion are exact. We 
have investigated quadratures squeezing, principal squeezing, quasiprobability distribution 
functions, quadrature distribution, phase distribution and phase variances. Generally, we 
have shown that the light obtained from the system exhibits Poissonian statistics and pro- 
vides squeezing in the framework of quadratures and principal squeezing. Moreover, it has 
been shown that when the values of A increase the period of the energy exchange between 
waveguides decreases. This fact leads to many interesting effects, such as the quadratures 
squeezing and the phase variances can exhibit collapse-revival and collapse-revival-subrevival 
phenomena, respectively. The generation of YSCS in the system has been analytically and 
numerically demonstrated and confirmed in all studied quantities. YSCS can be generated 
even if the linear interaction between the waveguides is neglected. Moreover, we have shown 
that the generation of these states most probably occurs when initially one of the modes is 
in the coherent state while the other is in vacuum state, provided that t\ = Ntt and N is 
a fraction of integer. Furthermore, the system can generate the statistical-mixture coherent 
states in dependence on the values of the interaction parameters. The nonclassical effects 
have been remarked in the behaviour of P(x,t). Finally, P(x,t) cannot include a complete 
information on the interference in phase space. 

Appendix 

In this appendix we give the derivation of the quasiprobability distribution given by (l2"Tj) . 
It is worth reminding that t\ = (m' + 1/2)tc, i.e. z = —1. We start by providing the 
following array: 

m = 1, 2,3,4,5,6, 7,8,..- 
l = -L( ni -l) = 0,0,1,3,6,10,15,21,- •• 

z l = 1,1,-1,-1,1,1,-1,-1,--- (41) 

From the information shown in ( 14T1) we can express the summation associated with the index 
rii in (1251) (for this case) as 

W(P, t X , s)=J2 {-l) ni z~^ n2 ~ l) [h{2 ni - n 2 )F(2m, n 2 ) - h{2 ni + 1 - n 2 )F{2 ni + 1, n 2 )\, 

ni,ri2 

(42) 
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where 



-2 . a . 2 \ a?{t)aT\t) ( 2 



F( ni ,n 2 ) = — -exp |/3M 1 , w z ^-n^-m 

7T(1 — S) \1 — S J n 2- \1 — S J \1 — S 

'(43) 

/i( ni - n 2 ) = exp[e{z n2 ~ ni - 1)]. (44) 

It is evident that, for these specified values of tx, when n\ — n 2 is even h{n\ — n 2 ) = 1 
otherwise h{n\ — n 2 ) = exp(— 2e). Moreover, the summation related to n 2 in (1421) can be 
similarly expressed as that over n\ and hence (jUD takes the form 

w{p,t x ,s)= E (-i) (ni+n2) 

{ j P(2n 1) 2n 2 )-F(2n 1 + l,2n 2 + l) 

+ exp[-2e][F(2m, 2n 2 + 1) - F(2m + 1, 2n 2 )]}. 
Now we show how the first summation can be evaluated in a closed form: 

J2(-l) {ni+n2) F(2n u 2n 2 ) = i ^(z)K+« 2 )[l + (-l) ni ][l + (-1)" 2 ]FK, n 2 ). (46) 

rai,n2 ni,n2 



Substitute (|4"3"j) into (JIB]) and use the generating function for the Laguerre polynomials [35] 

as 



exp(-a;Jfe)(l + A;) 17 = k n K~ n ( x )- ( 47 ) 



Therefore, the summation related to the index i%\ in (j46j) can be easily carried out and we 
obtain 



E ("l) (ni+ " 2) F(2n 1 , 2n 2 ) = ^ £ 1 ^ J [1 + (-l) n2 ]{[l + ^#] n2 

•i,ri2 "2 

(48) 



x exp 



22 

Now the summation over n<i can be straightforwardly evaluated and we obtain the following 
closed form expression: 

E (-l) ( " 1+n2) F(2n 1; 2n 2 ) = {exp - za x (t)| 2 ) + exp + ^(t)| 2 ) 

ni ,n 2 ^ 

(49) 

+2exp [-^- s (\(3\ 2 - l«i(*)l 2 )] cos (£(/3a#) 

Similar procedures have to be performed to obtain the other terms in (T4"5]) and then ff271) is 
obtained. 
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